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Abstract. Within a multiphase transport (AMPT) model, where the initial conditions are 
obtained from the recently updated HIJING 2.0 model, the recent anisotropic flow and 
suppression data for charged hadrons in Pb+Pb collisions at the LHC energy of ^/snn = 2.76 
TeV are explored to constrain the properties of the partonic medium formed. In contrast 
to RHIC, the measured centrality dependence of charged hadron multiplicity dN^/drj at LHC 
provides severe constraint to the largely uncertain gluon shadowing parameter s g . We find final- 
state parton scatterings reduce considerably hadron yield at midrapidity and enforces a smaller 
s g to be consistent with dN^/dr/ data at LHC. With the parton shadowing so constrained, 
hadron production and flow over a wide transverse momenta range are investigated in AMPT. 
The model calculations for the elliptic and triangular flow are found to be in excellent agreement 
with the RHIC data, and predictions for the flow coefficients « n (pT,cent) at LHC are given. 
The magnitude and pattern of suppression of the hadrons in AMPT are found consistent with 
the measurements at RHIC. However, the suppression is distinctly overpredicted in Pb+Pb 
collisions at the LHC energy. Reduction of the QCD coupling constant a„ by ~ 30% in the 
higher temperature plasma formed at LHC reproduces the measured hadron suppression. 



1. Introduction 

Collisions of heavy nuclei at the Relativistic Heavy Ion Collider (RHIC) [U [2j [31 0] and recently 
at the Large Hadron Collider (LHC) JSJ [6] have created a matter consisting of deconfined but 
strongly coupled quarks and gluons (sQGP). Evidence of this is provided by the hydrodynamic 
model analysis of elliptic flow data that requires an extremely small viscosity to entropy 
density ratio [8] and from the observed suppression of high-p-j- hadrons [9j [10] in central 
collisions relative to both peripheral and nucleon-nucleon collision. The suppression has been 
well established as due to parton energy loss via medium induced gluon emission The lost 
energy, whose magnitude depends on the parton density, reappears as soft hadrons [T2| IT3]. 

As the QGP formed in central Pb+Pb collisions at ^/snn = 2.76 TeV is at a density of 2.4 
larger and probes parton distribution at a smaller momentum fraction x than at RHIC, analysis 
of the recent data for bulk hadron production [144 [T5] , anisotropic flow [TB] and high-p-p hadron 
suppression at LHC |17t [P5] may provide crucial insight into the nuclear medium properties of 
the hot and dense QCD matter. 

While perturbative QCD can address only hard scatterings, formation of sQGP and abundant 
soft particle production suggest a highly nonperturbative physics which is not yet well- 
established within QCD. Thus to explore medium effects on particle production from soft to the 
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Figure 1. dN^/dr] distribution at RHIC (left panel) and LHC (right panel) in 0-5% central 
collisions. The predictions from HIJING (dashed line), HIJING+ZPC (dashed-dotted line), 
AMPT (solid line) are compared to data (solid circles) from BRAHMS [25] and ALICE [T4] . 



hard scattering regime relies on systematic inclusion of various stages of dynamical evolution of 
matter. 

2. The AMPT Model 

A MultiPhase Transport (AMPT) model [19] which combines the initial particle distribution 
from HIJING model [20J with subsequent parton-parton elastic scatterings via the ZPC parton 
cascade model and final hadron transport via ART allows a systematic study of hadron 
production. To investigate bulk and hard particle production, the AMPT model was modified to 
include the updated HIJING 2.0 version. In the two-component HIJING model [20] for hadron 
production, nucleon-nucleon collision with transverse momentum px transfer larger than a cut- 
off po leads to jet production. Soft interactions withpr < Po is characterized by an effective cross 
section <7 so f t . In HIJING 2.0 [2T] the Duke-Owens parametrization [22J of the parton distribution 
functions has been updated with the modern Gliick-Reya-Vogt (GRV) parametrization |23j . 
Since the gluon distribution at small momentum fraction x is much larger in GRV, instead of a 
fixed value for (po, <7 so f t ) as in HIJING 1.0, an energy dependent cut-off for po(y/s) and <7 so ft(\/s) 
are used to fit experimental data on total and inelastic cross sections in p + p/p collisions [21J. 

For the nuclear parton distribution function (PDF), HIJING employs the functional form 
f A (x, Q 2 ) = ARa(x, Q 2 )fa(x, Q 2 )-, where is the PDF in a nucleon. The nuclear modification 
factor of quarks and gluons (a = q,g) in HIJING 2.0 parametrization are [21] 

Rf(x, b) = h{x) - s q (b) (A 1 / 3 - I) ' 6 (1 - 3.5x - 5 ) exp(-x 2 /0.01), 

Rf(x,b) = /i(x)-s 9 (6)(A 1 / 3 -l) a6 (l-1.5x a35 )exp(-x 2 /0.004). (1) 

Here h(x) = 1 + 1.2 log 1 ^ 6 ^ (x 3 — 1.2x 2 + 0.21x), the impact parameter dependence of shadowing 
is taken as s a (b) = (5s a /3)(l — 6 2 /i? 2 ), and Ra ~ A 1 / 3 is the nuclear size. s q = 0.1 is fixed from 
deep inelastic scattering data and s g shall be estimated from fits to the measured charged particle 
pseudorapidity density in heavy ion collisions. We have used here the string melting version of 
the AMPT where parton recombination is employed for hadronization. At both RHIC and LHC 
energies, we consider the strong coupling constant a s = 0.33 and screening mass [i = 3.226 fm -1 
|24j that correspond to parton-parton elastic scattering cross section of a ~ 97ra 2 /(2 / u 2 ) ~ 1.5 
mb. 

3. Results and discussions 

3.1. Particle yield distribution 

Figured] shows the pseudorapidity distribution of charged hadrons, dN^/drj, for central collisions 
in Au+Au at ^/W^ = 200 GeV and Pb+Pb at y/sjfif = 2.76 TeV. The AMPT results 




Figure 2. dN^/dr] at mid-rapidity per participant nucleon pair as a function of (iVp ar t). 
The results are from AMPT calculations (triangles) with gluon shadowing s 9 =0.10-0.17 at 
RHIC (left panel) and with s g = 0.16 — 0.17 at LHC (right panel) as compared with the 
data (circles) from BRAHMS [25] and PHENIX [26J at RHIC and ALICE [T2] at LHC. 

are with gluon shadowing parameter of s g = 0.15 (at RHIC) and s g = 0.17 (at LHC) 
that agrees well with the measured dN^/drj distribution from BRAHMS at RHIC, and the 
dN ch /dr/ (\rj\ < 0.5) = 1601 ± 60 from ALICE at the LHC. The HIJING 2.0, without any 
final state interactions, predicts dN c h/drj (\r/\ < 0.5) = 706 ± 5 and 1775 ± 3 at RHIC and LHC, 
respectively. Energy dissipation and redistribution via partonic scatterings in subsequent parton 
cascade (i.e. in HIJING plus ZPC) result in the reduction of charged particle yield by ~ 15% at 
both RHIC and LHC inspite of ~ 2.4 times increase in the partonic density at LHC [27]. Finally 
scatterings in the relatively dilute hadronic matter (i.e. AMPT) lead to only a small decrease 
of particle abundances. 

Figure [2] shows the charged particle pseudorapidity density per participant pair, 
(dN c h/ drj) I ^ (iVp art ) /2) , as a function of centrality of collision characterized by average number 
of participating nucleons (A par t). The AMPT calculations are performed with a range of gluon 
shadowing parameter of s g = 0.10 — 0.17 at the RHIC energy and with s g = 0.16 — 0.17 at 
the LHC. This choice of s g leads to centrality dependence of charged particle multiplicity that 
agrees well within the experimental uncertainty of the BRAHMS [25] and PHENIX [26] data at 
RHIC. Due to abundant jet and minijet production in Pb+Pb collisions at LHC, the ALICE 
multiplicity results are quite sensitive to nuclear distortions at small x and provide a much 
stringent constraint on the gluon shadowing of s g ~ 0.17. We note that the estimated values 
of s g in AMPT are consistently smaller than the HIJING 2.0 [21] estimates of s g = 0.17 - 0.22 
(RHIC) and s g = 0.20 — 0.23 (LHC) which indicate the importance of final state interactions in 
precise determination of the nuclear shadowing of partons. 

3.2. Anisotropic Flow 

Partonic interactions within the initial spatially asymmetric collision zone drives uneven pressure 
gradients that manifests in anisotropic emission of particles leading to collective flow. The 
magnitude of this flow is characterized [281 [29] by the Fourier coefficients v n = (cos(n[</> — ^ n ])) 
as elliptic (^2), triangular (^3) and quadrangular (7J4) flow estimated with respect to their 
participant event planes Characterization can also be made via the correlations between the 
/c-particles with azimuthal angles <p\, . . . , <pk as v{ni, . . . , n^} = (cos(ni</>i + • • • + nk<ftk)), where 
77-1, . . . , nk are integers and the average is taken over all charged particles and events. This allows 
to construct the 2- and 4-particle correlations, v n {2} = \/v{n, —n} and v n {4}, respectively. 
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Figure 3. Top panels: px dependence of V2 (left panel), v% (middle panel) at various centralities, 
and centrality dependence of v n (n = 2, 3, 4) (right panel) at RHIC. The solid symbols are STAR 
[5U] and PHENIX [31] data. Bottom panels: pr dependence of v n at centralities of 0-5% (left 
panel), 30-40% (middle panel), and centrality dependence of v n (right panel) at LHC. The solid 
symbols are ALICE data j!6j. 



In Fig. [3] (top panel), the transverse momentum dependence of elliptic flow V2{2}{px) for 
charged hadrons at midrapidity is displayed at the RHIC energy at various centralities. It 
is seen that over a large px range the AMPT model is consistent with the STAR data |30j . 
Spatial fluctuations of participating nucleons lead to v% of similar magnitude as V2 for most 
central collisions in AMPT. While the triangularity in the collision geometry, and thus v^, is 
nearly insensitive to increase in collision peripherality (impact parameter), the Vi increases as 
the initial spatial eccentricity becomes more pronounced. We find AMPT also reproduces the 
measured W3{EP}(pt) (top-middle panel) as well as U4{EP}(pt) (not shown here) in the event- 
plane method at various centralities with the same cross section of a ~ 1.5 mb and shadowing 
s g = 0.15 at the RHIC energy as used/estimated above. The AMPT model predictions for the 
centrality dependence of v n {2} (n = 2,3,4) with momenta 0.2 < pr < 8 GeV/c are presented 
in Fig. [3] (top-right panel) at the RHIC energy. With collision centrality, the pT-averaged V2 
shows substantial variation, the changes are small in V3 and negligible in v±. 

For central Pb+Pb collisions at the LHC energy, the pt dependence of v n {2} as seen in Fig. [3] 
(bottom panel) exhibits similar magnitude and pattern as at RHIC energy. The enhanced parton 
density at the LHC energy that is expected to produce larger V2 appears to be compensated 
by faster expansion of the evolving medium due to larger flow. The centrality dependence of 
v n {2} (n = 2,3,4) is found to agree quite well with the ALICE data [TB]. As expected the 
smaller nonflow effects result in successively reduced elliptic flow in the event-plane f2{EP} and 
4-particle correlation U2{4} methods. 
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Figure 4. Raa for charged hadrons as a function of pt in central and peripheral collisions 
for Au+Au at RHIC (left panel) and Pb+Pb at LHC (right panel). The data are from 
STAR [TO] at RHIC and from ALICE [T7] and CMS [TO] at LHC. The AMPT results are 
with strong coupling constant a s = 0.33 at RHIC and LHC and with a s = 0.24 for central 
collisions at LHC 



3.3. High-pT charged hadron suppression 

The study of bulk hadron production in conjunction with that for hadron spectra provide crucial 
information on the parton-medium interactions where high-p-j- partons suffer energy loss which 
in turn produce soft hadrons. The suppression of hadrons at high pt due to medium effects in 
heavy ion collisions is quantified by the nuclear modification factor 

d 2 N AA /d V dp T 

Raa(pt) = (N con ) ffiim/dndpr (2) 

which refers to the ratio of particle yield in heavy ions (A + A) to that in p + p reference spectra, 
scaled by the total number of binary nucleon-nucleon (NN) collisions (N co u) = (Taa)o~^ ■ In 
absence of initial and final state nuclear medium effects Raa(pt) = 1 by construction. The 
nuclear thickness function (Taa) and the inelastic NN cross section are calculated within 
the HIJING 2.0 that employs Glauber model for the distribution of initial nucleons with a 
Woods-Saxon nuclear density. 

The nuclear modification factor Raa for charged hadrons is shown in Fig. UJ For central 
collisions at RHIC and LHC, Raa(pt) < 1 which suggests appreciable suppression of particles 
relative to NN reference. With the nuclear shadowing parameter s g = 0.15 constrained from 
dN c h/dr] data in Au+Au collisions, the AMPT results describe the magnitude and pattern of 
the RHIC suppression data [TO]. The success of AMPT at RHIC suggests that the initial state 
shadowing, the final state scattering and the parton energy loss is consistent with the formation 
and evolution of the medium at the RHIC energy. 

In peripheral Pb+Pb collisions at LHC, the Raa for is nearly constant at about 0.7 over 
a large pt range. Here, the QGP even if formed, should have a small volume and short lifetime. 
In central Pb+Pb collisions at LHC, the rise and fall pattern exhibited by Raa up to pr ~ 6 
GeV/c (in data and model) is similar to that at RHIC. The subsequent rise of Raa at LHC is 
mainly due to the harder unquenched pQCD jet spectra. However, in contrast to ALICE and 
CMS data, the AMPT calculations show even more pronounced suppression at pt > 2 GeV/c. 



This suggests that the LHC medium with a factor of 2.4 increase in density over RHIC is, in 
fact, more transparent [27} [32]. 

It may be noted that at both RHIC and LHC, the Raa was evaluated in the AMPT with 
the same value for the QCD coupling constant a s = 0.33 and screening mass \i = 3.226 fm _1 . 
As the screening mass depends on temperature as p, = gT = \J\^a s T [33], the parton-parton 
elastic scattering cross section used in AMPT becomes a ~ 97ra 2 /(2^ 2 ) r* 9a s /(8T 2 ). Using the 
scaling relation between the initial entropy density Sj and the final particle multiplicity [34(. 155] 
S{ ~ (7.85/Tj^4j_) dN c h/dy, and considering a QGP of massless gas of light quarks and antiquarks 
so that Si ~ 4ej/(3Tj) (with energy density a « (21/30)vr 2 T. 4 ), allows one to estimate the initial 
temperature T and thereby the parton scattering cross section a from the measured particle 
yield. In — 5% central Au+Au and Pb+Pb collisions at ^/snn = 0.2 and 2.76 TeV, the 
measured dN ch /dy « 687 [26] and 1601 [33] result in T rj 320 and 436 MeV respectively, at a 
proper time Tj = 1 fm/c. With the above choice of a s = 0.33, the estimated a « 9a s / (8T 2 ) « 1.4 
mb at RHIC turns out to be similar to the value employed in AMPT that reproduces the RHIC 
suppression data shown in Fig. 03 On the other hand, the higher Tj at LHC leads to a ~ 0.76 
mb. Alternatively, if fi is set constant at 3.226 fm from RHIC to LHC, this smaller a then 
gives a s ~ 0.24 at LHC. With this reduced a s , we find from Fig. [l]the AMPT results for Raa 
(open circles) in central Pb+Pb collisions at LHC is in good agreement with the suppression 
data. This is a clear indication of thermal suppression of the QCD coupling constant at the 
higher temperature reached at LHC. 

4. Conclusions 

In summary, within the AMPT model that is updated with the HIJING 2.0 initial conditions 
for parton distributions, we study the properties of the medium formed in heavy ion collisions 
at RHIC and LHC energies. We find energy redistribution via parton scatterings reduces the 
final hadron multiplicity and thereby enforces a smaller gluon shadowing for the hadron yield 
to be consistent with the data. The magnitude and trend of the flow coefficients, t>2, i>3, t>4 as 
a function of pr and at various centralities are in excellent agreement with the measurements 
at RHIC and LHC. The parton energy loss in AMPT is found to describe the charged hadron 
suppression over a large px range at RHIC. However, using the same QCD couplings a s and 
parton scattering cross section from RHIC to LHC, jet quenching in AMPT is distinctively 
overpredicted relative to the measurements at central Pb+Pb collision at LHC. A reduction of 
a s by ~ 30%, consistent with the higher temperature of the plasma formed at LHC, agrees with 
the measured suppression. 
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